Observation of vortex dipoles in an oblate Bose-Einstein condensate 
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We report experimental observations and numerical simulations of the formation, dynamics, and lifetimes 
of single and multiply charged quantized vortex dipoles in highly oblate dilute-gas Bose-Einstein condensates 
(BECs). We nucleate pairs of vortices of opposite charge (vortex dipoles) by forcing superfluid flow around a 
repulsive gaussian obstacle within the BEC. By controlling the flow velocity we determine the critical velocity 
for the nucleation of a single vortex dipole, with excellent agreement between experimental and numerical 
results. We present measurements of vortex dipole dynamics, finding that the vortex cores of opposite charge 
can exist for many seconds and that annihilation is inhibited in our highly oblate trap geometry. For sufficiently 
rapid flow velocities we find that clusters of like-charge vortices aggregate into long-lived dipolar flow structures. 

PACS numbers: 03.75.Kk, 03.75.Lm, 67.85.De 



Vortex dipoles consist of a bound pair of vortices of op- 
posite circulation and may exist in both classical and quan- 
tum fluids. Although single vortices carry angular momen- 
tum, vortex dipoles can be considered as basic topological 
structures that carry linear momentum 111 in stratified or two- 
dimensional fluids. Vortex dipoles are widespread in classi- 
cal fluid flows, appearing for example in ocean currents yfl 
and soap films yfl, and have been described as the pri- 
mary vortex structures in two-dimensional chaotic flows [1]. 
In superfluids, the roles of quantized vortex dipoles appear 
less well established. Given the prevalence of vortices and 
antivortices in superfluid turbulence Jj-la], the Berezinskii- 
Kosterlitz-Thouless (BKT) transition [7[|, and phase transition 
dynamics ||- IT], a quantitative study of vortex dipoles will 
contribute to a broader and deeper understanding of superfluid 
phenomena. The realization of vortex dipoles in dilute Bose- 
Einstein condensates (BECs) is especially significant as BECs 
provide a clean testi ng grou nd for the microscopic physics of 
superfluid vortices I12H1411 . In this paper we present an ex- 
perimental and numerical study of the formation, dynamics, 
and lifetimes of single and multiply charged vortex dipoles in 
highly oblate BECs. 

Numerical simulations based on the Gross-Pitaevskii equa- 
tion (GPE) have shown that vortex dipoles are nucleated when 
a superfluid moves past an impurity faster than a critical veloc- 
ity, above which vortex shedding induces a drag force [ 15 , Kj . 
Vortex shedding is therefore believed to be a mechanism for 
the breakdown of superfluidity IU71I18I1 . Experimental studies 
of periodic stirring of a BEC with a laser beam have measured 
a critical velocit y fo r the onset of heating and a drag force on 
superfluid flow lll9l EoTl . and vortex phase singularities have 



been observed in the wake of a moving laser beam 12 U 122 1. 



deterministic ally nucleated by causing a highly oblate, har- 
monically trapped BEC to move past a repulsive obstacle. We 
measure a critical velocity for vortex dipole shedding, and 
find good agreement with numerical simulations and earlier 
theory |f23[l. Experimentally, the nucleation process exhibits 
a high degree of coherence and stability, allowing us to map 
out the orbital dynamics of a vortex dipole. We find that vor- 
tex dipoles can survive for many seconds in the BEC without 
self-annihilation. We also provide evidence for the formation 
of multiply charged vortex dipoles. 

The creation of BECs in our lab is described elsewhere 



However, a microscopic picture of vortex dipole formation 
and the ensuing dynamics has not been established experi- 
mentally. In the work reported here, single vortex dipoles are 



| ITj, |24|, |25|] . In the experiments reported here, we begin with 
a BEC of 2 x 10 6 atoms in a highly oblate harmonic trap. 
Our axially symmetric trap is created by combining a red- 
detuned laser light-sheet trapping potential with a magnetic 
trapping potential, producing a BEC with an 11:1 aspect ra- 
tio and a Thomas-Fermi radius of 52 pun radially, as shown 
in Fig.[Tta,b). The BECs are additionally penetrated by a fo- 
cused blue-detuned laser beam that serves as a repulsive ob- 
stacle; the beam has a Gaussian 1/e 2 radius of 10 /an and is 
initially located 20 /jm to the left of the minimum of the har- 
monic trap as shown in Fig. [TJc). To nucleate vortices we 
translate the harmonic potential in the horizontal (x) direction 
at a constant velocity until the obstacle ends up 14 fim to the 
right of the harmonic trap minimum. At the same time, the 
height of the obstacle is linearly ramped to zero as shown in 
Fig-[TJd,e), allowing us to gently create a vortex dipole that is 
unaffected by the presence of an obstacle or by heating due to 
moving the obstacle through the edges of the BEC where the 
local speed of sound is small. For our conditions, the healing 
length at trap center is ~ 0.3 pirn, which is approximately the 
size of vortex cores in our trapped BECs. After a subsequent 
variable hold time t/, we remove the trapping potential and ex- 
pand the BEC for imaging, causing the vortex cores to expand 
relative to the condensate radius such that they are optically 
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FIG. 1. Illustration of the BEC initial state and the experimental se- 
quence, (a) Side-view phase-contrast image and (b) axial absorption 
image of a BEC in the highly oblate harmonic trap in the absence of 
the obstacle. Lighter shades indicate higher column densities inte- 
grated along the line of sight. Our axial and radial trapping frequen- 
cies are oj- = 2n x (90 Hz) and a> r = 2n x (8 Hz), respectively, (c) 
BEC initial state with the obstacle located at x s = -20 yum relative 
to the BEC center. (d,e) The maximum repulsive potential energy 
of the obstacle is Uq ~ 1.2/io (where /io ~ 8fiw ; is the BEC chem- 
ical potential) and is held constant during evaporative cooling. It is 
ramped down linearly as the trap translates; relative to the trap cen- 
ter, the beam moves from position x, = -20 yum to x s = 14 fjm over 
a variable sweep time t s . The BEC is then held in the harmonic trap 
for a variable time % prior to expansion and absorption imaging. 



resolvable. An example axial absorption image is shown in 
the left-most inset image of Fig. [2] 

In Fig. [2] we plot the average number of vortices observed 
as a function of the trap translation velocity v s . In our ex- 
perimental procedure we find a ~ 30% likelihood of a single 
vortex occurring during condensate formation even prior to 
trap translation; see Ref. Ill ill for further details. This gives 
a non-zero probability of observing a single vortex for the 
lowest translation velocities in Fig. [2] even when flow with- 
out drag is expected. The results of zero-temperature and 
finite-temperature c-field numerical simulations 12711 are also 
shown — for simulation details see Ref. 82511 . There is good 
agreement between simulation and experimental results, and 
we identify a critical velocity v c for vortex dipole formation 
between 170 /iin/s and 190 yimi/s for N c -2x 10 6 atoms and 
temperature T — 52 nK. Recently, Crescimanno et al. 12311 
have calculated the critical velocity for vortex dipole forma- 
tion in a 2D BEC in the Thomas-Fermi regime. By using the 
nonlinearity and chemical potential of our 3D system in their 
2D expression, we estimate a critical velocity of 200 yum/s. 
For our conditions, the maximum speed of sound at the trap 
center is calculated to be c ~ 1700 /iin/s; our measurements 
show that v c ~ 0.1c, consistent with previous measurements 
of a critical velocity for the onset of a drag force If20ll . 

In an axi-symmetric trap such as ours, a vortex dipole co- 
incides with a meta-stable state of superfluid flow with poten- 
tially long lifetimes 1 28 -3(1. The vortices exhibit periodic or- 
bital motion, a 2D analogue of the dynamics of a single vortex 
ring l22ll3lll . To observe these dynamics we nucleate a single 
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FIG. 2. Average number of vortex cores observed for a range of 
translation velocities v s with N c ~ 2 x 10 6 atoms in the BEC. Ex- 
perimental data points (black diamonds) are averages of 10 runs 
each, with error bars showing the standard deviation of the obser- 
vations. Numerical data for N c = 2 X 10 6 at a system temperature of 
T = 52 nK, corresponding to the experimental conditions, are indi- 
cated by triangles joined with dotted lines. Fewer atoms and lower 
temperatures reduce the critical velocity; as an example, open circles 
show the results of numerical simulations using N c = 1.3 x 10 6 and 
T = 0. For v s < 170 yum/s, remnant vortices arising spontaneously 
during condensate formation (see Ref. fTltl ) account for the baseline 
of ~ 0.3. For v s ~ 190 yum/s, the inset images show a single pair 
of vortices having formed in the experiment (left) and simulations 
(right). Because the vortex core diameters are well below our imag- 
ing resolution, the BEC is expanded prior to imaging. Similarly, the 
vortex cores in the unexpanded numerical image are barely visible; a 
10-yum-wide inset provides a magnified scale for the core size. Above 
v s ~ 200 /jm/s, multiply-charged vortex dipoles are observed. The 
critical velocity calculated using the methods of Ref. l23ll is indicated 
by the vertical dashed line. 



vortex dipole and hold the BEC for variable time tf, prior to 
expansion and imaging, with results shown in Fig. [3ja). The 
repeatability and coherence of the vortex nucleation process 
is clear: in back-to-back images with increasing %, the vortex 
positions and orbital dynamics can be followed and the dipolar 
nature of the superfluid flow is microscopically determined. 
These measurements also determine the direction of super- 
fluid circulation about the vortex cores, analogous to the case 
of single vortices 13211 : the image sequence shows counter- 
clockwise fluid circulation about the vortex core in the upper 
half of the BEC and clockwise circulation in the lower half. 
The orbital dynamics were also examined in zero-temperature 
GPE simulations, as shown in Fig. [3jb) |25ll . and the exper- 
imental and numerical data are in good agreement regarding 
vortex dipole trajectories, as shown in Fig.[3jc). The lifetime 
of a single vortex dipole is much longer than the first orbital 
period of ~1.2 s (see below), although after the first orbit the 
vortex trajectories become less repeatable from shot-to-shot. 
However, it is the large-scale flow pattern of the first orbit 
that is perhaps most indicative of the qualitative connection 
between 2D superfluid and classical dipolar fluid flows. 
For trap translation velocities well above v c we observe 
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FIG. 3. Sequences of images showing the first orbit of vortex dipole dynamics, (a) Back-to-back expansion images I2al from the experiment 
with 200 ms of successive hold time between the 1 80-/jm-square images. This data sequence was taken using an obstacle height 3.15 times 
larger than that used for the data of Fig.[2] as we find this gives the highest degree of repeatability and the least sensitivity to beam displacement, 
(b) 62-yum-square images from numerical data obtained for conditions similar to the data of sequence (a), but for a temperature of T = 0. The 
apparent vortex core size is smaller in the simulations as we show in-trap data rather than expanded BECs. (c) Black and dark gray small 
circles show average positions x, and y v of each of the two vortices from 5 sequences of experimental data, each sequence using a procedure 
identical to that of sequence (a). The larger circle around each average position point represents the standard deviation of the vortex positions 
at that specific hold time, and is calculated from the 5 images obtained at that time step. A continuous dipole trajectory from sequence (b) 
is re-scaled to the Thomas-Fermi radius of the expanded experimental images, and superimposed as solid lines on the experimental data; no 
further adjustments are made for this comparison, (d) Similar to the sequence of (b), but for a faster translation velocity with which a doubly 
charged vortex dipole is formed. Images are spaced in time by 120 ms as the vortices orbit more quickly, (e) An experimental image in which 
a doubly charged vortex is formed; see text for further explanation. (f,g) Similar to (d,e) but for a triply charged dipole and 100 ms between 
images in (f). 



the nucleation of multiply-charged vortex dipoles both experi- 
mentally and numerically, as shown in Fig.[3{d-g). Viewed on 
a coarse scale the ensuing dynamics are consistent with that of 
a dipole comprised of a highly charged vortex and antivortex. 
On a fine scale, particularly in numerical data, we see loose 
aggregations of singly quantized vortices with the same circu- 
lation at the two loci of vorticity in the dipolar flow. In the ex- 
perimental images obtained at higher sweep speeds, many in- 
dividual vortices are often not clearly resolvable for the short 
hold times shown. Nevertheless, the data resemble charac- 
teristics of highly charged dipoles and suggest the formation 
of many vortices because (i) the apparent vortex core sizes be- 
come larger, (ii) the orbital time period for the dipole structure 
is shorter, as expected for higher numbers of cores and faster 
superfluid flow, and (iii) multiple individual vortex cores are 
observable for longer hold times. Although we have not per- 
formed an exhaustive analysis of these states, we present these 
results to bring attention to these interesting metastable super- 
fluid vortex structures. 

While it is often assumed that in a finite-temperature envi- 
ronment, vortices of opposite circulation will attract and an- 
nihilate each other at close distances, this is not necessarily 
the case: vortices may approach each other so closely that 
they appear to coalesce — see for example the sixth image 
of Fig. [3{a) with f/, = 1 s — and yet still move away from 



each other after the encounter. In Fig. [4] we show measure- 
ments of the average number of vortices observed with vari- 
ous hold times after nucleating a vortex dipole, from which we 
conclude that singly and doubly quantized vortex dipoles may 
exhibit lifetimes of many seconds, much longer than a single 
orbital period. With such a strong trap asymmetry, the vortex 
lines are relatively impervious to bending (331 and tilting J34ll . 
and annihilation is suppressed because vortex crossings and 
reconnections are inhibited. 

Although vortex shedding is the microscopic mechanism 
for the breakdown of superfluid flow and the onset of a drag 
force, the nucleation of a vortex dipole does not imply im- 
mediate superfluid heating since dipoles are coherent struc- 
tures that are metastable and do not immediately decay into 
phonons. We estimate the maximum energy of a vortex dipole 
in the BEC is kg X 0.45 nK/atom for our system fl9tl . so that 
at T — 52 nK the temperature would increase by less than 
0.5% upon self-annihilation of a single dipole. This temper- 
ature change would be very difficult to measure, and micro- 
scopic observations of vortex dipole formation and dynamics 
are therefore complementary to heating and drag force obser- 
vations. 

In conclusion, we have demonstrated controlled, coher- 
ent nucleation of single vortex dipoles in oblate BECs. The 
dipole dynamics reveal the topological charges of the vortices 
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FIG. 4. Number of vortices remaining after dipole nucleation as a 
function of hold time t h , averaged over 17 realizations per data point. 
The circles show conditions for which a singly charged dipole is cre- 
ated, while the squares show data from a faster sweep where four 
cores (a doubly charged dipole) preferably occur. The error bars rep- 
resent statistical uncertainty (as in Fig. [2} rather than counting un- 
certainty; however, for the doubly charged dipole case there is addi- 
tional uncertainty in counting vortex cores because we do not always 
resolve 4 well-defined cores at the earlier hold times. The vortex 
cores can clearly persist for times much longer than the dipole or- 
bital period of ~1.2 s for a singly quantized dipole, and ~ 0.8 s for 
a doubly quantized dipole. The inset figures show experimental im- 
ages of singly quantized vortices seen after 4 seconds of hold time in 
the harmonic trap. 

and show that the dipole is a long-lived excitation of super- 
fluid flow. Sufficiently rapid translation of the harmonic trap 
causes vortices with identical charge to aggregate into highly 
charged dipolar vortex structures that exhibit orbital dynamics 
and long lifetimes analogous to singly charged vortex dipoles. 
This suggests that both dipole structures and macro-vortex 
states are readily accessible in highly oblate and effectively 
two-dimensional superfluids. 
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AUXILIARY MATERIAL FOR: Observation of vortex dipoles in an oblate Bose-Einstein 
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Here we provide auxiliary details regarding our experimental methods and theoretical techniques for studying 
and modeling vortex dipoles in Bose-Einstein condensates, as reported in the main paper [T. W. Neely et al.]. 



Experimental methods: In previous work fU, Q], we de- 
scribed our methods to create 87 Rb BECs in the F = 1, nif = 
-1 state in a time-averaged orbiting potential (TOP) trap, 
which was combined with optical dipole potentials from a 
blue-detuned laser. To study BECs in a highly oblate ge- 
ometry, we further compress the trap in the vertical (z) di- 
mension by adding a horizontally propagating, 1-W, sheet of 
red-detuned laser light of wavelength 1090 nm. The resulting 
trap has radial and axial trapping frequencies of (u> r , u> z ) = 
2jt x (8,90) Hz, with the radial frequency dominantly de- 
termined by the TOP trap portion of the potential. Radio- 
frequency evaporative cooling in this highly oblate harmonic 
trap achieves condensation at a critical temperature of T c ~ 90 
nK, and produces a BEC of N c = 2 + 0.5 x 10 6 atoms at a fi- 
nal temperature of T ~ 52 nK. The BEC chemical potential 
in the harmonic trap is ~ %hcj z , placing our BECs well 
within the three-dimensional (3D) regime. Radial and axial 
BEC images are shown in Fig. l(a,b) of the main paper. 

The repulsive gaussian obstacle used to nucleate vortex 
dipoles is formed by a blue-detuned 660-nm laser beam propa- 
gating along the axial (z) direction of the trap. This beam has a 
Gaussian 1/e 2 radius of 10 /an at the BEC, much smaller than 
the 52 pm Thomas-Fermi radius of the BEC, and the height of 
the repulsive potential has an energy of approximately 1.2^q. 
The beam is turned on prior to condensation, and is aligned 
20 + 3 /mi to the left of the trap center as shown in Fig. 1(c) of 
the main paper. After the BEC forms in this annular trap, the 
harmonic potential is translated in the horizontal (x) direction 
at a constant velocity until the beam position ends up 14 jjm to 
the right of the trap center and the translation then stops. Trap 
translation is accomplished by changing a magnetic bias field, 
a process that we find to be more repeatable and stable than 
moving the beam with piezo-electric transducers. To initiate 
trap translation, the position of the trap center is smoothly ac- 
celerated over ~5 ms before reaching a constant translation 
velocity for the duration of the translation stage. This brief 
acceleration stage ensures that residual "sloshing" of the BEC 
due to excitation of center-of-mass oscillations is minimized; 
experimentally we find the amplitude of such oscillatory mo- 
tion to be less than 1 fim. Coincident with trap displacement, 
the beam intensity is ramped to zero as shown in Fig. l(d,e) 



of the main paper. After a variable hold time f/,, we enhance 
vortex visibility by expanding the BEC radially for 15 ms in 
the presence of the laser sheet with the TOP trap off, followed 
by free expansion for an additional 41 ms. An example axial 
absorption image is shown in the left inset image of Fig. 2 of 
the main paper. 

Simulations for T = 0: At zero temperature we simulate the 
experiment using a 3D Gross-Pitaevskii equation (GPE). We 
find an initial state for the combined harmonic-plus-obstacle 
trap with 2 x 10 6 atoms by integrating the GPE in imag- 
inary time. The potential due to the obstacle laser beam 
is U s (x,y,t) = U Q (t)e\p{-2([x-x (t)] 2 +y 2 )/wl), where 
i/o(0) = 93 ha) r = kg x 35 nK ~ /j.q, and wo = 10 yum. 
We assume that Us is constant over the axial position co- 
ordinate z, and x and y are radial position coordinates with 
(x, y, z) = (0, 0, 0) corresponding to the harmonic trap center. 
We model the vortex dipole nucleation procedure by translat- 
ing the obstacle while leaving the harmonic trap center fixed 
at (x,y,z) = (0,0,0). The obstacle laser beam moves along 
the x-axis from position xo(0) = -20yum to Xo(t s ) = 14/zm 
during time t s and at a constant velocity v s = (34 / um)/f J . As 
in the experiment, the laser intensity, and thus Uo(t), linearly 
ramps to zero over time t s . We then continue the simulations 
to observe the T = dynamics of the condensate and any vor- 
tices that were nucleated. We vary v s over the range 133 to 
256 yum/ s in steps of 8 /an/s. 

We emphasize that the simulations described above do not 
directly correspond to the experimental sequence, where the 
harmonic trap is translated rather than the obstacle. While 
the appropriate timescale for strictly adiabatic dynamics is 
<d~ 1 = 20ms «: t a , a shorter acceleration stage is experimen- 
tally found to be suitable for the trap motion used, as men- 
tioned above. Thus our simulations are equivalent to applying 
an instantaneous Gallilean boost to the condensate at t — so 
that the condensate is stationary in the translating frame. We 
refer to this as the "adiabatic" regime of trap acceleration. 

We have performed additional simulations to investigate the 
"impulse" regime where the harmonic trap is instantly accel- 
erated to the sweep velocity. In this case we perform the sim- 
ulation in the frame co-moving with the trap — the trans- 
formation to this frame results in an initial phase gradient 
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across the BEC causing it to initially move with the obsta- 
cle. This procedure excites a center-of-mass oscillation, typ- 
ically of amplitude 10 fim or more near the critical velocity, 
much larger than observed in the experiment. These numeri- 
cal observations thus suggest that the experiment is far closer 
to the "adiabatic" regime results presented in the main paper. 
We note that simulations in the impulse regime also show a 
significantly lower critical velocity (~ 80 yum/s), and multi- 
ple vortex pairs are much more readily nucleated at lower trap 
translation velocities. This appears to be due to relatively high 
velocities between the obstacle and the oscillatory motion of 
the BEC center-of-mass (due to instantaneous trap accelera- 
tion); factors other than the constant velocity of the trap must 
therefore be considered. A complete model of the experiment 
would need to incorporate details of the trap acceleration that 
were not obtained. However, given the excellent agreement 
between the adiabatic regime of simulations and the full set 
of experimental results, including the experimental observa- 
tions showing minimal BEC sloshing after trap acceleration, 
we conclude that our simulation methods are appropriate for 
comparisons with the experimental data. 

Simulations for T > 0: While the experiments are per- 
formed at the lowest attainable temperatures, there is still a 
non-negligible thermal fraction present. For this reason we 
have performed simulations of the vortex nucleation process 
using c-field simulation techniques J3{]. Broadly this is an ex- 
tension of the Gross-Pitaevskii equation to finite temperatures 
where the lowest energy modes of the system have a large 
thermal occupation and it is reasonable to treat them classi- 
cally. 



Our procedure is as follows. We find an initial finite- 
temperature state for the system by evolving the stochastic 
Gross-Pitaevskii equation yl |4[] to thermal equilibrium for an 
input chemical potential of kg x 37 nK corresponding to a 
condensate number of 2.0 x 10 6 atoms, and a temperature of 
52 nK as measured in the experiment. Once we have gen- 
erated a number of initial states for a given temperature, we 
apply the same simulation procedure as described above for 
the zero-temperature GPE. We perform our simulations in the 
harmonic oscillator basis, and implement a basis cutoff at 130 
nai r . This is necessary in order to accurately account for the 
modification of the condensate wave function by the obstacle 
beam — in doing so we include a rather large range of ex- 
cited states that are not highly occupied in the classical field 
ifr, hence violating the usual validity condition for classical 
field simulations. However, this is unavoidable, and we ex- 
pect that it will not have significant effects on the nucleation 
of the vortex dipoles, although it may affect their subsequent 
dynamics. 
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